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We study the dilute fermion gas with pairing between two species and unequal concentrations in a 
harmonic trap using the mean field theory and the local density approximation. We found that the 
system can exhibit a superfiuid shell structure sandwiched by the normal fermions. This superfluid 
shell structure occurs if the mass ratio is larger then certain critical value which increases from the 
weak-coupling BCS region to the strong-coupling BEC side. In the strong coupling EEC regime, the 
radii of superfluid phase are less sensitive to the mass ratios and are similar to the case of pairing 
with equal masses. However, the lighter leftover fermions are easier to mix with the superfluid core 
than the heavier ones. A partially polarized superfluid can be found if the majority fermions are 
lighter, whereas phase separation is still found if they are heavier. 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.90.-|-q 



I. INTRODUCTION 

Recently experimental progress has raised strong in- 
terest in studyingthe superfluid pairing in the ultra- 
cold Fermi gases Through the Fcshbach resonance 
Q , the effective interaction between atoms can be tuned 
over a wide range. This technique led us to investigate 
the crossover between the condensation of weak-coupling 
Cooper pairs and the Bose-Einstein condensation (BEC) 
of strong-couphng tightly bound pairs Q. More recently, 
two experimental groups extend a step further by con- 
trolling the polarization of the two states Fermi gases Q 
and open a new era for studying the imbalanced Fermi 
gases, which is a topic related to many interesting areas 
from condensed-matter physics, nuclear physics, to quark 
matter f^. 

Stimulated by these excellent experiments of the im- 
balanced Fermi gases, there are intense theoretical works 
in the past two years. For a homogeneous system, the 
smooth crossover is known to be destroyed when the pop- 
ulations of the fermions are unequal [q, 0, H, @] ■ In an 
inhomogeneous case, e.g. in a harmonic trap potential, 
phase separation occurs near resonance with the super- 
fluid is at the trap center surrounding by the normal 
phase [13, mi, [H, El, [3, ■ Finite temperature phase 
diagrams and density profiles are also studied in this im- 
balanced system recently ^1, [13, [H, [H, [13l ■ 

Feshbach resonances between different atomic species 
have also been reported [IT'l and the pairing of unequal 
masses fermions has received theoretical interest. Focus- 
ing on the pairing of ^Li-^'^K mixtures, the phase diagram 
was studied both in a homogeneous system [23| and in 
a trap potential [23j . Other works on a general unequal 
mass ratios, including the analytic results at unitary limit 
[2^ and the evolution of the tricritical point [1^ have also 
been studied recently. 

In this paper, we study the density profiles for unequal 
mass fermions in harmonic traps at zero temperature. 



Within the local density approximation, or Thomas- 
Fermi approximation, we solve the BCS gap equation 
self-consistently at fixed total number of particles and 
polarization. We analyze the density profiles from un- 
polarized to highly imbalanced phases. Due to the large 
number of possible atoms (e.g. ^H, ^He* [1^, ^Li, ^"K, 
^''Sr [27,], and ^''^Yb [11]), we study general mass ratios 
from 0.1 to 10 and particular {kpa)^^ = —0.5, 0, and 1.0 
as the representative scenarios for weak-couping BCS, on 
resonance, and strongly paired BEC regimes. We found, 
on the weak coupling BCS side, the superfluid can form a 
shell structure sandwiched between normal fermions for 
large mass ratios. At resonance, the normal fluids can 
be either partially or completely polarized, depending on 
the mass ratios. On the BEC side, the radn of superfluid 
core are insensitive to the mass ratios. A partially po- 
larized superfluid can be found if the majority fermions 
are lighter, whereas phase separation is still found if they 
are heavier. We also report the axial density profiles of 
the population difference and discuss the significant dif- 
ferences compared to the pairing with equal masses. 

Mismatch Fermi surfaces may result a different ground 
state than phase separations, in particular the so-called 
Fulde-Ferrel-Larkin-Ovchinnikov (FFLO) phase [13] • In 
this paper, we leave out this possible FFLO phase for 
future investigation [soj . 

The remainder of this paper is organized as follows: In 
Sec. II we briefly review the mean-fleld approximation 
for the dilute two states of fermion atoms with unequal 
masses. We present our numerical results along with dis- 
cussion, in Sec. III. Finally, in Sec. IV, we conclude with 
a briefly summary. 



II. FORMALISM: 

We start from the two-component fermion system 
across a wide Feshbach resonance which may be described 
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by an effective one-cliannel Hamiltonian: 

k,o- 

+ 5 51 4+q,h.4'-q,;Ck',iCk,/. , (1) 
k,k',q 

wliere ^o-lk) — h^hF' /2m„ — g the bare coupling 
strength, and the index a runs over the two species [h 
and I). Within the BCS mean field approximation, the 
excitation spectrum in a homogeneous system for each 
species is (see e.g. [3l| for details) 
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where —h = l. In order to include the pairing between 
unequal masses, we use rrih {mi) to represent the mass of 
each component with heavier (lighter) mass and rewrite 
the dispersion relation of Eq. ([2]) as 



ft^fc^ 7 - 1 
4to,- 7+1 



h 



(3) 



with the reduced mass nir, the mass ratio 7 — mh/mi, 
the chemical potential difference h = (fih — M;)/2 and the 
average chemical potential fx = {^h + /^i)/2. 

For the fermion gases in a harmonic trap, the system 
can be treated as homogeneous locally if the number of 
particles are sufficiently large. Within this local density 
approximation, the densities at position f depends only 
on the local effective chemical potentials 



(4) 



where Vh.i{r) is the trap potential for the species h or 
I. To simplify our presentation we shall consider an 
isotropic trap where the potential depends only on the 
radius r. 

The density profiles can be found via 
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with / is the Fermi function. The total number of parti- 
cles N — J d^rNs{r) and the polarization of the system 
is defined as 
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Note that the polarization is positive (negative) for a 
system with the majority are heavier (lighter). 

Now the pairing field A depends on position also. In 
the local density approximation, it obeys an equation 
similar to the homogeneous case |3lj : 

(Pk 
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For a given s-wave scattering length a, we solve equations 
(O, (O and ^ self-consistently for fixed total number 
of particles N and polarization P. The solutions of the 
"gap equation" ([5]) may not be unique and the physi- 
cal solution is determined by the condition of minimum 
free energy among the multiple solutions. The detailed 
procedure can be found in reference [lo| . 

To avoid extra complications, we confine ourselves to 
the case where the trapping potential is harmonic. We 
further assume that it is identical for the two species. 
Thus V{r) — ^ah,ir'^ with ah = cti = a. This would 
occur, for example, if both species are trapped magnet- 
ically and if the fermions have identical magnetic mo- 
ments, or when the fermions are trapped optically with 
different lasers of appropriate intensities and detunings. 
In this case, within the local density approximation, the 
density profiles for the same polarization P but differ- 
ent total particle numbers are related to each other via 
simple scaling. This allow us to present our results in a 
manner which is independent of the total particle num- 
bers (see below). When the trap potentials are unequal, 
many other scenarios can occur. We shall leave the in- 
vestigation of this more complicated case to the future. 



III. RESULTS 

In this section, we investigate the density profiles for 
various polarizations, mass ratios and coupling strengths 
from positive detuning BCS superfluid to negative de- 
tuning EEC side. With the aid of density profiles, we 
plot the superfluid phase diagram in a harmonic trap for 
different mass ratios. Finally, we also evaluate the axial 
density profiles of the population difference and discuss 
the effect due to the mass ratios at different coupling 
strengths. 

We begin by examining the radial density profiles for 
the mass ratios 7 = 2.0, and 5.0 with various polariza- 
tions from 0.8 (majority are heavier) to -0.8 (majority 
are lighter). The densities N{r) are normalized to the 
total density at the trap center and distances are normal- 
ized to the Thomas-Fermi radius of the non-interacting 
majority gas of the same atom number. In Figs. [T]-[31 
we plot the radial density profiles at three different cou- 
pling strengths with {kpa)~^ — (the unitary limit), 1.0 
(strong coupling BEC), and -0.5 (weak-coupling BCS). 

In the unitary limit, the system exhibits a superfluid 
core surrounded by normal fermions for 7 = 2.0 [Fig. 
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[Ija)-(e)]. The density profiles are qualitatively the same 
between positive and negative polarizations. For exam- 
ple, for the large polarizations P = ±0.8 [Fig. [^a) and 
(e)], we found that the superfluid core of the majority are 
lighter is only slightly larger than that of the majority are 
heavier. However, for P > there are small amounts of 
minority in the normal region but it is completely po- 
larized for the case of P < 0. The above results are in 
accordance with Fig. 2 of reference At the inter- 

face between the superfluid and the normal fluid, the two 
phases are in equilibrium and the local polarization of the 
latter phase should be given by the same value obtained 
from our earlier bulk calculations there. The asymmetry 
between the positive and negative polarizations becomes 
more obvious for larger mass ratio, e.g., 7 = 5 in Fig. 
[Uf)-(j). In Fig. mg) and (h), the superfluid is sand- 
wiched between two regions with normal fermions. The 
inner core is a normal state with heavy particles being the 
majority, whereas the outer normal shell consists purely 
the lighter species The superfluid shell extends to- 
ward the trap center as the polarization decreases from 
large positive value. At equal populations [Fig. [Hh)], we 
found the system is still phase separated into normal po- 
larized fermions and unpolarized superfluid regions. This 
is quite different compare to the equal or small mass ratio 
{e.g., 7 = 2.0 here) cases which the system is superfluid 
for the entire trap [33| . 

The existence of this normal fluid core can be under- 
stood as due to the fact that this normal fluid has a 
higher number density than the paired superfluid at this 
mass ratio [13], and hence it "sinks" towards the center 
of the trap where the trap potential is the minimum (re- 
call that we have assumed ah — on). We note also here 
that, as shown in Fig [l] the total number density is a 
decreasing function of r. In fact, the results of [24] then 
indicate that the critical mass ratio above which for the 
occurrence of a normal state core with a superfluid shell 
is 7* = 3.9 at resonance (see Fig. 3 of reference [24j). 

In Fig. [21 we plot the radial density profiles for 7 = 2.0 
and 5.0 with (kpa)^^ = 1.0. At this coupling strength, 
the radii of the superfluid core are found to be only very 
weakly dependent on the mass ratio or the sign of polar- 
ization. However, the phases are very different depend- 
ing on whether the majority particles are the heavier or 
lighter ones. We found that all particles are paired for a 
unpolarized case (P = 0). When the majority fermions 
are heavier (P > 0), the system contains an unpolar- 
ized superfluid core with a surrounding outside shell of 
normal fluid consisting of the majority particles alone. 
Thus phase separation occurs in this case between the 
completely paired superfluid and the completely polar- 
ized normal phase. On the other hand, when the major- 
ity is light (P < 0), the equally paired superfluid core 
at the trap center is surrounded by a partially polar- 
ized superfluid mixture phase, with a normal shell oc- 
cupying the outermost region of the trap. This normal 
shell is again found to be completely polarized. In this 
strongly paired EEC regime, all minority particles exist 



only within the bound pairs and thus within the super- 
fluid phase(s) only. Note also that the couphng strength 
is now strong enough so that the superfluid always has 
a higher density than the normal phase, and the super- 
fluid shell structure found at resonance no longer appears 
here. 

In the BCS regime, since the pairing interaction is 
weak, the density profiles are most easily understood by 
first considering the normal state density profiles. We 
note here that due to the different mass ratios and un- 
equal populations, these two profiles have very different 
dependence on the position r. In the presence of the weak 
interaction, pairing would occur only at locations where 
the local populations arc almost equal, or when the den- 
sity is sufficiently large so that the effective interaction 
is sufficiently strong. This results in a rich structures 
in the density profiles as we shall show. In Fig. [3l we 
plot the radial density profiles for the coupling strength 
(/cpa)"^ = —0.5 and mass ratios 7 = 2.0 and 5.0. The 
symmetry between the positive and negative polariza- 
tion is lost even for 7 = 2.0 [Fig. [3Ka)-(e)]. At large 
populations imbalances {e.g., P = ±0.8 in Fig. [3]), a su- 
perfluid core still forms for majority are lighter but only 
normal fermions exist in the trap for majority are heav- 
ier. Similar to the unitary limit, the system exhibits a 
superfluid core surrounded by normal fermions for a wide 
range of polarization at 7 = 2.0. But at larger 7 [Fig. 
|3l[f)-(j)], the heavier normal fermions prefer to stay in 
the trap center and the superfluid phase is outside this 
normal phase. We also found that the superfluid shell 
grows as the population of lighter fermions increase [Fig. 
[Hg)~(i)]- This superfluid shell structure finally reaches 
the trap center and all normal fermions are forced into 
the outside shell. We would like to emphasize that the 
shell structure of the normal fermions are not the same 
as the cases in the unitary limit. In Fig. [^g) and (i), 
the outside shell of the normal fermions are fully polar- 
ized but both components of fermions exist here for this 
weak-coupling side. 

In Fig. ^h), we found that the density profile con- 
tains two regions of superfluid and polarized fermions. 
As mentioned, a stable superfluid occurs when either 
the density of each component of fermion is large or the 
size of Fermi surfaces of these two fermion states are ap- 
proximately equal. In Fig. 21 we enlarge the Fig. [3Kb) 
and plot the density profiles of both components of the 
non-interacting fermions also. It shows that the super- 
fiuid shell just occurs near the crossing of these two non- 
interacting density profiles since the superfluid gap is eas- 
ier to open near this point [2^. On the other hand, the 
density is large enough such that the superfluid phase is 
again stable at the trap center. 

In view of the rich structure on the BCS side, we turn 
next to the phase diagram of the superfluid in a trap 
potential in this regime. We plot the phase diagrams in 
Fig. [5] for several different mass ratios with the same 
coupling strength {kpa)~-^ — —0.5. In Fig. [SJa) (7 = 1), 
the superfluid phase extends to the entire trap at the un- 
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polarized system (P = 0) and is symmetric about this 
point. For 7 > 1, the fermion paired between differ- 
ent masses, the maximum radius of the superfluid phase 
moves toward negative polarization where the majority 
of the system are lighter. The symmetric phase diagram 
does not hold anymore. At increasing mass ratio (7 > 3), 
we found that the superfluid at the trap center is not sta- 
ble for the positive polarization. Instead, the superfluid 
phase occurs at some finite radius for the parameters we 
studied in Fig. [ntb)-(e). The shell structure of the su- 
perfluid is present for a wide range of polarizations as 
the mass ratio 7 increases. In Fig. [Hb), 7 = 2.0, the 
superfluid phase is split into two regions by the normal 
fermions and vice versa for < P < 0.4, correspond- 
ing to the concentric structure we discuss in Fig. Eljb) 
above. Note that, the mass ratios in Fig. [H^b), (c), and 
(d) correspond to the mixtures of ■^He*-^Li, ^H-^Li and 
6Li-40K. 

In Fig. [HI we also plot the superfluid phase for three 
different coupling strengths with 7 — 5.0. The superfluid 
shell structure still exists at the unitary limit [Fig. ^h)] 
but only for positive polarizations. On the BEG side 
[Fig. EIc)], the symmetry of the superfluid phase about 
the equal population almost regains. However, there is 
a region where the superfluid and leftover fermions are 
mixed with each other for P < 0. For P > 0, there 
is only unpolarized superfluid phase and the system is 
phase separated. 

Lastly we would like to examine the axial density pro- 
files of the population difference defined as : 

Nd{z) = J dxdy[Nh{rO - Ni{r)] . (9) 

In Fig. [71 we plot these results for (kpo)^^ — —0.5, 0.0 
and 1.0 and the mass ratios (7' = rrimajoity /rrirmnority) 
between majority and minority species from 0.1 to 10.0 
at polarization |P| = 0.2. 

On the BCS side [Fig. [7{a)] and the unitary limit 
[Fig. [7{b)], the horizontal segments of each curve Nd{z) 
correspond to the unpolarized superfluid. For the cases 
of majority lighter (7' < 0.2), we found that the minority 
(with heavier mass) crowd into the trap center and have 
higher population than the majority locally. Thus, the 
axial density proflles show the dip structure in the trap 
center. In the other limit, 7' > 5.0 (the majority are 
heavier), the peak of the axial density profile appears 
right at the trap center. 

At the unitary hmit [Fig. EUb)], the axial density pro- 



files are similar except for 7' = 5.0 and 10.0. This is 
the results where the heavier fermions have larger pop- 
ulations near the trap center and the superfluid shows a 
shell structure as we discuss in Fig. [T] 

On the strong coupling EEC side [Fig. [TJc)] , the struc- 
tures of all of the curves are similar for different mass 
ratios. The radii of the superfluid core in this coupling 
strength are not sensitive to 7' either. The only differ- 
ence is that Nd{z) are smooth at the phase boundary for 
the cases of the majority are lighter but there exist kinks 
for the cases of the majority arc heavier . It is consistent 
with the results in Fig. [2] that the lighter fermions are 
easier to mix with superfluid than the heavier ones. 



IV. CONCLUSION 

We have studied the two-component fermion system 
with unequal masses and populations across Feshbach 
resonance. For the pairing with unequal species, the su- 
perfluid can be sandwiched by the normal fermions and 
form a superfluid shell structure. This superfluid shell 
structure is easier to observe on the weak-coupling BCS 
side or at the unitary limit with large mass ratios (i.e., 
in the systems of ^°K~-^'^^Yh, ^Li-^OK, or ^Li-S^Sr^). 
For a given mass ratios, this superfluid shell extends to- 
ward the trap center as the population of lighter fermions 
increase. 

In the strong coupling BEG regime, the superfluid 
phase is less sensitive to the mass ratios and is similar to 
the case of pairing with equal masses. However, lighter 
leftover fermions are easier to penetrate into the super- 
fluid core than the heavier ones. At coupling strength 
(fci?a)^^ = 1.0, phase separation between the completely 
polarized superfluid and the completely polarized normal 
phase occurs for the majority are heavier. On the other 
hand, the equal paired superfluid core at the trap center 
is surrounded by a partially polarized superfluid mixture 
phase and a normal shell occupying the outmost region 
of the trap. 
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FIG. 1: (Color online) The radial density profiles for heavy 
fermions (solid lines) and light fermions (dashed lines) at the 
unitary limit {(kpa)'^ = 0). The dotted lines are the total 
density profiles. The mass ratios 7 are equal 2.0 for plots (a)- 
(e) and 5.0 for plots (f)-(j). The vertical scale is normalized 
to the total density of fermions at trap center (A^s(O)). ttf 
is the Thomas-Fermi radius of the normal majority cloud. 
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FIG. 3: (Color online) The radial density profiles on the BCS 
side with {kpa)'^ = —0.5. The parameters are the same as 
in Fig. 1 
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FIG. 4: (Color online) The radial density profiles for heavy 
fermions (solid lines) and light fermions (dashed lines) with 
parameters (kpa)'^ — —0.5, 7 = 2.0, and P = 0.2. The dot- 
ted lines are extrapolated density profiles of each component 
as if they were non-interacting. Solid circles represent the dif- 
ference between the two components of fermions. Note that, 
Nd is slightly less than zero near the edge of the trap. 
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FIG. 5: (Color online) Superfluid phase diagram for 
{kFa)~^ = —0.5 in the trap. The shaded area represents 
the superfluid phase. 
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FIG. 6: (Color online) Superfluid phase diagram for 
{kpay^ = -0.5 (a), 0.0 (b) and 1.0 (c). The mass ratio 
7 = 5.0. The vertical shaded area shows unpolarized super- 
fluid phase and the horizontal shaded area (in plot (c)) shows 
partially polarized superfluid phase. 
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FIG. 7: (Color online) Axial density profiles for {kpa) ^ = 
-0.5 (a), 0.0 (b) and 1.0 (c) with |P| = 0.2. Solid lines 
correspond to the majority heavier (7' > 1) and dashed lines 
corresponding to the majority lighter (7' < 1). 



